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ABSTRACT 

We present the results from a Chandra pilot study of 12 massive galaxy mergers selected from Galaxy 
Zoo. The sample includes major mergers down to a host galaxy mass of 10 11 Mq that already have 
optical AGN signatures in at least one of the progenitors. We find that the coincidences of optically 
selected active nuclei with mildly obscured (Nh )$ 1-1 x 10 22 cm -2 ) X-ray nuclei are relatively common 
(8/12), but the detections are too faint (< 40 counts per nucleus; / 2-10 keV < l-2xl0 -13 ergs -1 cm -2 ) 
to reliably separate starburst and nuclear activity as the origin of the X-ray emission. Only one merger 
is found to have confirmed binary X-ray nuclei, though the X-ray emission from its southern nucleus 
could be due solely to star formation. Thus, the occurrences of binary AGN in these mergers are rare 
(0-8%), unless most merger-induced active nuclei are very heavily obscured or Compton thick. 

Subject headings: galaxies: active — X-rays: galaxies 


1. INTRODUCTION 

Major mergers are a key component of current mod- 
els for galaxy formation in a ACDM Universe. Merg- 
ers can disrupt the star-forming gas and stellar disks of 
the progenitors, trigger a powerful burst of star forma- 
tion, and reshape the remaining stellar content into a 
bulge. Perhaps with a small time delay, the supermas- 
sive black holes may feed on gas from the destabilized 
or destroyed disk, injecting energy in the form of radia- 
tion or kinetic outflows that sweep the remnant clear of 
dust and gas. First proposed by Sanders et al. (1988), 
this picture directly links the triggering of active galac- 
tic nuclei (AGNs) phases to the destructive potential in- 
duced by a merger. Recent semi-analytic models and 
hydrodynamic simulations have adopted this scenario to 
explain the fueling of AGNs and the red spheroidal rem- 
nants that are difficult to reproduce without some kind of 
“AGN feedback” (Springel et al. 2005; Di Matteo et al. 
2005; Hopkins et al. 2006, 2008; Somerville et al. 2008). 
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In principle, major mergers carry with them two black 
holes, both of which may be accreting and be visible 
as distinct AGN during a phase of abundant gas avail- 
ability that a major, gas-rich merger represents. Yet the 
evidence associating AGN phases with major mergers re- 
mains contested (De Robertis et al. 1998; Malkan et al. 
1998; Schmitt 2001; Pierce et al. 2007; Georgakakis et al. 
2009; Gabor et al. 2009; Schawinski et al. 2011). Large 
optical surveys using Sloan Digital Sky Survey (SDSS) 
data have found ~3.6% of spectroscopically confirmed 
AGNs are in closed binaries (~5-100 kpc separation; 
Liu et al. 2011). The DEEP2 survey also found that 
binary AGN exist in ~2.2% (2/91) of red galaxies 
with type 2 Seyfert optical spectra (Gerke et al. 2007; 
Comerford et al. 2009). However, optical surveys can 
easily miss obscured AGNs especially in merger systems 
where the gas is driven toward the center through dis- 
sipation (e.g., Hopkins et al. 2008). X-ray surveys are 
needed to identify the more highly obscured systems 
(Nh > 10 20 cm -2 ). We know of only a small number 
of binary AGN resolved directly using X-ray observa- 
tions (e.g., Komossa et al. 2003; Guainazzi et al. 2005; 
Hudson et al. 2006; Bianchi et al. 2008; Foreman et al. 
2009; Comerford et al. 2011; Fabbiano et al. 2011). The 
intrinsic frequency of binary AGN phases has not been 
observationally constrained, as the separation of individ- 
ual X-ray sources is not possible in high-redshift sources 
and there has been no systematic search for such sys- 
tems in known mergers. A study of the host galaxies of 
185 nearby (z < 0.05) BAT AGNs by Koss et al. (2011) 
found that these hard X-ray selected AGNs are preferen- 
tially found in massive galaxies with large bulge-to-disk 
ratios and large supermassive black holes. This may im- 
ply that the frequency of binary AGNs is higher in mas- 
sive mergers. 

In order to quantify the intrinsic frequency of double 
AGNs in the local universe, we embarked on a study of 
the presence of binary AGNs and their dependence on the 
mass of the host galaxies. The results from the present 
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survey represent a pilot effort as the sample is comprised 
of only the most massive galaxies in the Galaxy Zoo 
merger sample. Thus, this paper aims to quantify the 
intrinsic frequency of double AGNs in the mass limit 
down to ^10 11 M 0 using a study of 12 merging galaxies 
with the Chandra X-ray observatory. Throughout this 
paper, we adopt H 0 = 71 km s _1 Mpc -1 , VLm = 0.3, and 
n A = o.7. 

2. SAMPLE SELECTION 

The initial parent sample was created from 3003 merg- 
ers identified via visual inspection by citizen scientists 
taking part in the Galaxy Zoo project (Lintott et al. 
2008; Darg et al. 2010a, b). To date, this is the largest un- 
biased sample of mergers in the local Universe (0.005 < 
z < 0.1). From this catalog, we whittled the sample 
down to only 328 major mergers (i.e. , with a mass ra- 
tio of 3:1 or less) in which the primary galaxy has a 
mass 13 > 10 11 Mq; we also required that the galaxies 
have SDSS spectroscopic data with signal-to-noise ratio 
greater than 3 and that at least one of the nuclei shows an 
AGN signature based on [O III]/H/? and [N II] /Ha nar- 
row line ratios (Baldwin et al. 1981; Kewley et al. 2006). 
In addition, the AGN must be relatively luminous, with 
-^[O hi] > 10 41 erg s" 1 (yellow points in Figure 1). This 
last cut is to ensure that the AGN will be luminous 
enough for detection in the X-rays. Most of the objects 
in this final sample are LINERs or AGN and star-forming 
composites, from which we selected 12 that spanned the 
full starburst-composite-LINER/AGN range. The pri- 
mary nuclei of the selected sample (yellow dots circled 
in red) also have roughly even distributions in the com- 
posite but not extreme starbursts (40%) and the LINER 
(53%) regimes of Figure 1 and these mergers have pro- 
jected nuclear separations between 5 to 14 kpc. The 
optical line ratios are taken from Oh et al. (2011). The 
12 galaxies in our sample are listed in Table 1; for conve- 
nience, we will identify the targets as Galaxy Zoo (GZ) 
objects throughout this paper. For clarity, we will refer 
to the merging galaxies as either “mergers” or “galax- 
ies” and the individual progenitors of these mergers as 
“nuclei” for the remainder of this paper. Figure 2 is a 
collage of the SDSS images of the sample. All of these 
mergers show disturbed morphology indicative of tidal 
disruptions. 

3. OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

The 12 mergers were observed with Chandra between 
November 2010 and May 2011 (PI: Schawinski). Each 
pair of galaxies was observed in a single exposure of 
4.9 ks. For each merger, the more massive primary nu- 
cleus was placed at the nominal aim point of the ACIS-I3 
chip. The progenitors of the mergers are close enough that 
the secondary nuclei were also within the same chip. 

The data were reduced using CIAO version 4.3 and 
CALDB version 4.4.3. The data reduction followed pro- 
cedures outlined in the Science Analysis Threads for 

13 The stellar masses of the galaxies were calculated following 
the methodology outlined in Schawinski et al. (2010). Briefly, mea- 
surements from the five SDSS photometric bands were fitted to a 
library of model star formation histories generated from Maraston 
(1998, 2005) stellar models. Stellar masses are measured by finding 
the minimum of the x 2 statistic in the parameter space probed. 


ACIS imaging data on the CIAO web page 14 . For each 
of the 12 detected nuclei, we extracted counts in the soft 
(0.5-2 keV) and hard (2-8 keV) bands. Two of these 
nuclei were detected having only two counts in the soft 
band, but none in the hard band. Therefore, only 10 
nuclei have valid hardness ratios (HRs 15 ; Table 1) for es- 
timating spectral properties. For the rest of this paper, 
only the 10 nuclei with valid HRs will be discussed. 

An HR analysis was performed using the latest version 
of FTOOLS released as part of HEASoft 6.11. Due to 
the low number of counts measured from our sample, we 
followed Teng et al. (2005) and used the measured HR 
to estimate a photon index (r) by assuming a redshifted 
power law model modified only by Galactic absorption. 
The nominal HR and the estimated photon indices are 
tabulated in Table 1. The X-ray fluxes were then esti- 
mated using PIMMS by assuming the HR-derived photon 
indices and the count rates from the observations. 

4. DISCUSSION 

The shapes of the X-ray spectra differ for AGNs, ob- 
scured AGNs, starbursts, and AGN-star forming com- 
posites. Typically, unobscured AGNs have spectra that 
are well-represented by a power law with photon index 
of ~1.8. Obscuration affects the lower energy (<2 keV) 
photons more readily than the higher energy photons and 
thus flatten or harden the AGN spectra. Starburst spec- 
tra are dominated by emission in the lower energies, but 
low-mass X-ray binaries tend to have relatively fiat spec- 
tra. Composite objects generally have softened spectra 
compared to simple AGN spectra due to the soft-energy 
contribution of the starburst. 

Given that these mergers contain optically selected 
AGNs, it is unsurprising that eight of the 10 nuclei have 
HRs that are consistent with the canonical spectral shape 
of unobscured AGNs (T ~ 1.7 — 2.1). As many as five 
could be steeper (GZ IE, GZ 4S, GZ 5S, GZ 7SW, and 
GZ 9S), as if star formation is a significant contribu- 
tor, though the errors in HR and T allow for unobscured 
AGN values. Similarly, three nuclei (GZ 2N, GZ 9N, 
and GZ 11S) have nominally flat spectra, implying dom- 
inance from star formation or obscured nuclear activity. 
However, the errors in HR cannot rule out unobscured 
AGN as the source of the X-ray emission. Finally, two 
nuclei (GZ 3S and GZ 10S) have flat or inverted pho- 
ton indices (T < 1.45 after accounting for the measure- 
ment errors) , suggesting some level of obscuration. If we 
assume a power law with T fixed at 1.8, the observed 
HRs imply column densities (Nh) < 10 21-22 cm -2 (Ta- 
ble 1). At these column densities, the HR estimates of 
the 2-10 keV luminosity are reliable to within ~40% 
(Teng & Veilleux 2010). These columns do not suggest 
the presence of Compton-thick nuclei though there re- 
mains a possibility of leaky, Conrpton-thick absorbers. 

4.1. Starburst Contamination 

Of the 12 mergers in the sample, one has no X-ray 
detection (GZ 8) and only one (GZ 9) exhibits binary 
X-ray nuclei (Figure 3). The remaining 10 mergers have 
one detected nucleus each. In GZ 5, GZ 7, and GZ 12, 

14 http:/ /cxc. harvard.edu/ciao/threads/index.html 

15 HR = / | / . where H and S are the total counts in the hard 
and soft bands, respectively. 
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the X-ray-detected nucleus is not the one with an op- 
tical AGN classification, so in that sense they are dou- 
ble nuclei. In addition, the detected southern nucleus of 
GZ 10 has extended soft X-ray emission (Figure 3), sug- 
gesting a contribution from star formation. This raises 
the question of whether more of the detections might be 
contaminated by star formation. 

To explore this possibility, we compare star formation 
rates derived from the SDSS u band luminosities follow- 
ing Hopkins et al. (2003) with those derived from the 2- 
10 keV luminosity following Ranalli et al. (2003) in Fig- 
ure 4. When compared with the SDSS M-band derived 
star formation rates (Figure 4), the X-ray derived star 
formation rates of four nuclei (GZ IE, GZ 5S, GZ 7SW, 
and GZ 9S) have unconstrained lower limits. While the 
nominal X-ray derived star formation rates are above the 
line of equality implying the presence of AGNs, we can- 
not rule out the possibility that the X-ray emission can 
be accounted for solely by star formation in these four nu- 
clei. The X-ray luminosities of the remaining nuclei are 
above those expected from star formation even after the 
consideration of the <40% error in the calculation of the 
X-ray luminosity, consistent with additional contribution 
to the X-ray luminosity by nuclear activity. Accounting 
for the error bars, the southern nucleus in GZ 9 may also 
be dominated by star formation, suggesting GZ 9 does 
not contain an AGN pair. 

4 . 2 . Compton-thick Nuclei 

Three nuclei of the 12 SDSS-selected mergers are not 
detected in X-rays (GZ 5N, GZ 7NE, and GZ 12E). We 
already know these have optically identified AGN com- 
ponents, so it is unclear whether the non-detections are 
due to faint AGNs (two of the three have the highest 
redshifts in our sample) or Compton-thick AGNs. If 
we assume these are faint AGNs, a power law model 
with r = 1.8 and mild absorption from the Milky Way 
places upper limits to the luminosity of these objects. 
Not accounting for intrinsic absorption, the 2-10 keV lu- 
minosity for GZ 5N is < 5 x 10 39 erg s _1 cm -2 and 
< 5 x 10 4 ° erg s" 1 cm" 2 for GZ 7NE and GZ 12E. In 
the case of the Compton-thick AGNs, the optical signa- 
ture is coming from the much larger scale narrow- and 
broad-line regions while the X-ray is sensitive to the small 
scale emission from the black hole itself. The presence 
of undetected obscured nuclei would affect our statistics 
of the frequency of binary AGNs. It is unlikely that 
all of the secondary nuclei contain Compton-thick X-ray 
sources, unless an obscured phase is common to mergers 
(unlike isolated AGNs). Even without a merger-induced 
obscured phase, the number of heavily obscured AGNs 
is comparable to the number of less obscured AGNs 
(Treister et al. 2009) ; the presence of Compton-thick nu- 
clei remains a possibility. 

While the individual detected nuclei have too few 
counts for spectral fitting to definitively establish 
whether Compton-thick AGNs are present, we consid- 
ered the cumulative rest-frame photon distribution of the 
detected nuclei in the hard band. We compared this ob- 
served distribution with the expected photon distribu- 
tions from unobscured AGNs and from Compton-thick 
AGNs. In the former case, we assumed a single unab- 
sorbed power law; in the latter case, we assumed a power 
law with an iron emission line at 6.4 keV with an equiv- 


alent width of 1 keV, a typical signature of Compton- 
thick AGNs. For both cases, the total photon counts 
were normalized to be the same as the total detected 
counts. In Figure 5 we plot the cumulative distribution 
of the detected photons in our sample. There is no clear 
distinction between the observed distribution with either 
model. In fact, the two-tailed Kolmogorov-Smirnov (K- 
S) test statistics for the two cases are nearly identical. 
As a sanity check, we compared the two modeled distri- 
butions with each other and there is a clear difference 
at the ~80% confidence level. Therefore, we cannot rule 
out the possibility that Compton-thick AGNs are present 
at the level that we are able to detected these sources. 

4 . 3 . U/LIRGs in Formation? 

In theoretical models of galaxy mergers (e.g., 
Hopkins et al. 2008), luminous and ultraluminous in- 
frared galaxies (U/LIRGs) represent a stage that merg- 
ers go through before the formation of elliptical galax- 
ies. Initially, tidal torques enhance star formation and 
black hole accretion. Then in the final coalescence of the 
galaxies, massive inflows of gas trigger starbursts with 
strengths similar to those inferred for U/LIRGs. 

The mergers in our sample appear to be the predeces- 
sors to U/LIRGs in this evolutionary picture. The X-ray 
luminosities estimated for our mergers are approximately 
10 times lower than those observed in most U/LIRGs, 
but are consistent with the lower end of the range mea- 
sured in LIRGs (Teng & Veilleux 2010; Lehmer et al. 
2010; Iwasawa et al. 2011). This implies mergers in our 
sample are in the earliest stages of interaction, where the 
growth of the central black hole has not yet peaked. 

The incidence of binary AGNs in U/LIRGs is also 
rare. The Revised Bright Galaxy Survey (RBGS; 
Sanders et al. 2003) is a flux-limited sample of U /LIRGs 
from the IRAS All Sky Survey. Of the 629 extragalactic 
objects with 60 /im flux greater than 5.24 Jy, 86 are in- 
teracting galaxies that are visually similar to our sample 
in the optical (i.e. close binaries). Of these, 32 have high- 
quality X-ray data from either Chandra or XMM-Newton 
that is sensitive to the presence of an AGN. Not account- 
ing for the presence of undetected Compton-thick nuclei, 
only 3% (1/32) of the RBGS sources with X-ray data 
show binary X-ray nuclei (NGC 6240; Komossa et al. 
2003). This is consistent with the 0-8% (0-1 out of 12) 
we observe in our modest SDSS sample. 

4 . 4 . Frequency of Binary AGNs in SDSS Mergers 

From the very short snapshots of our study, we have 
found that coincidence of optically selected active nucleus 
with mildly obscured X-ray nucleus is relatively common 
(8/12). Given the faint detections, these snapshots are 
too short to place strong limits on the absence of AGN 
in the undetected galaxies, so it is difficult to comment 
on the frequency of binary active nuclei. However, we do 
detect a pair of X-ray nuclei in GZ 9, implying that this is 
uncommon unless the second nucleus is heavily obscured. 
In that instance, the most likely scenario would be that 
all nuclei are obscured. That is, either binary nuclei are 
uncommon, or merger nuclei in general have a high prob- 
ability of being heavily obscured. The latter possibility 
cannot be addressed by the current sample. To do bet- 
ter, we will need to increase the exposure times, expand 
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our merger sample for better statistics, and include a 
sample of major mergers for which there are no optically 
detected nuclei. Another natural follow-up would be to 
extend the study to a similarly selected sample with a 
lower mass limit to examine the dependence of binary 
AGNs on the mass of the host galaxies. 
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Fig. 1. — Emission line diagnostic diagram used for the selection of the Chandra Galaxy Zoo sample. The grey scale represents the 
complete Galaxy Zoo sample and the green points are the Galaxy Zoo AGN selected using narrow line diagnostics from Schawinski et al. 
(2010). The dashed curve shows the empirical separation between purely star- forming galaxies and the composite region of the diagram as 
determined by Kauffman et al. (2003, Ka03). The solid curve is the theoretical extreme starburst line of Kewley et al. (2001, keOl) beyond 
which the dominant source of ionization must be due to something other than star formation. The straight line demarcates the empirical 
AGN-LINER separation in Schawinski et al. (2007, S07). The yellow points are the Galaxy Zoo mergers that meet the criteria of mergers 
having mass ratios of at least 3:1 and having one of the progenitors with a mass > 10 11 Mq with significant emission line detections; 
most are LINERs. The blue points are the secondary nuclei with SDSS optical spectra. The Chandra observed nuclei (both primary and 
secondary) are circled in red with their GZ identifier from Table 1 labeled. Our Chandra sample covers the full range of [O III]/H/3 and 
[N II] /Hq emission line ratios for composite and LINER objects. They are also representative of the merger distribution in the composite 
and LINER areas of this Baldwin-Phillips-Terlevich (BPT) diagram (Baldwin et al. 1981). 
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Fig. 2. — Composite gri images of the Chandra - observed sample from SDSS Data Release 7. Each frame is labeled with the Galaxy Zoo 
identification listed in Table 1 and the horizontal bar represents angular distance of 20 arcseconds. 
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Fig. 3. — False color X-ray images of two interesting objects in our sample, GZ 9 (left) and GZ 10 (right). The red represent the 0.5-2 keV, 
green the 2-6 keV, and blue the 6-8 keV emission. The contours are from SDSS i band images. The raw X-ray images were smoothed with 
a 0'/5 Gaussian, the width of the nominal point spread function of Chandra. The binary nuclei in GZ 9 are both detected in the X-ray, 
though the southern nucleus is dominated by soft X-ray emission. In GZ 10, the X-ray emission shows east-west extension which may be 
due to star formation in addition to an obscured AGN. The flatness of the X-ray spectrum implies a column density ~ 5 x 10 21 cm - 2 if 
we assume a canonical power-law photon index of T ~ 1.8. 
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Fig. 4. — Comparison of star formation rates derived from the 2-10 keV luminosity (Ranalli et al. 2003) and SDSS u band luminosity 
density (Hopkins et al. 2003). Only the detected nuclei with hardness ratios, and thus estimated X-ray luminosities, are plotted. Errors 
are 1 — a. The arrows indicate poorly constrained negative error bar for three nuclei whose HR lower limits approach — 1, where the Nh 
and r values become degenerate for any value of HR. The solid line is the line of equal star formation rates. Each of the detected nuclei 
is labeled corresponding to the identification in Table 1. All of the detected nuclei have X-ray luminosity above that expected from star 
formation. 
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Fig. 5. — Cumulative rest- frame photon distribution in the hard band (2-8 keV) for the 12 detected nuclei. The black solid histogram 
represents the combined detected counts from all 12 detected sources in our sample. The blue dotted curve represents the expected 
cumulative distribution assuming an unobscured AGN model where a single unabsorbed power law represents the source of the AGN 
emission. The red dashed curve shows the distribution for a Compton-thick AGN model where the emission is represented by a power law 
plus a 6.4 keV iron emission line with an equivalent width of 1 keV. While there is a clear difference between the two modeled distributions, 
when compared with the observed distribution, neither can be shown to be the preferred model by a two-tailed K-S test. Therefore, we 
cannot rule out contribution from Compton- thick AGNs at present. 


